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■ The existing observational data on possible variations of fundamental physical constants (FPC) 
^ . confirm more or less confidently only a variability of the fine structure constant a in space and 

time. A model construction method is described, where variations of a and other FPCs (includ- 
es ■ ing the gravitational constant G) follow from the dynamics of extra space-time dimensions in 

■ the framework of curvature-nonlinear multidimensional theories of gravity. An advantage of this 
method is a unified approach to variations of different FPCs. A particular model explaining the 
observable variations of a in space and time has been constructed. It comprises a FRW cosmology 

^ with accelerated expansion, perturbed due to slightly inhomogeneous initial data. 



^3 ; 1 Introduction 



oo 

The problem of possible variations of the fundamental physical constants (FPC) in time and space 
\ is one of the most challenging problems of modern physics, directly related to the central problem of 
unification of all interations. It traces back to Dirac's and Eddington's famous papers of the 1930s 
O ■ and since then gains much attention in both theoretical and experimental studies. 

However, to date, a variability of only one FPC has been revealed by observations more or less 
confidently, it is the fine structure constant a . The analysis of absorption spectra of various ions in 
the radiation of distant quasars, performed in the recent years (above all, from the data obtained at 
the Keck telescope on the Hawayian islands), has led to a conclusion that a is changing with time, 
so that in the past it was slightly smaller than now (the relative change 5a /a is about 10~ 5 [1]). 
In 2010, an analysis of new data obtained at the VLT (Very Large Telescope), located in Chile, 
and their comparison with the Keck data led to a conclusion on spatial variations of a, i.e., on its 
dependence on the direction of observations. The VLT observations in the Southern part of the 
celestial sphere gave values of the parameter a in the past slightly larger than now. This anisotropy 
has a dipole nature [2,3] and has been termed "the Australian dipole" [4]. The dipole axis is located 
at a declination of —61 ±9° and at a right ascention of 17.3 ± 0, 6 hours. The deflection of a value 
at an arbitrary point r of space from its modern value a$ , measured on Earth, is 

5a/a = (1.10 ± 0.25) x 10" 6 r cos ip, (1) 

where ip is the angle between the direction of observation and the dipole axis, while the distance 
r is measured in billions of light years. The confidence level of this result (as compared with a 
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"monopole" model where values of a are the same in all directions) has been estimated as 4.1a. A 
more detailed discussion of the observational data can be found, e.g., in [3]. 

Let us also mention the laboratory experimental data on possible FPC variations in the modern 
epoch. The tightest constraints on a variations have been obtained by comparison of readings of 
atomic clocks using optical transitions in Al and Hg ions (without using cesium clocks that have 
become classic) [5]: (da/dt)/a = (—1.6 ± 2.3) x 10~ 17 per year. This result is of the same order of 
magnitude as the tightest constraints obtained previously from an isotopic composition analysis of 
the decay products in the natural nuclear reactor that operated in the Oklo region (Gabon) about 
2 billion years ago. Unlike the laboratory data, the Oklo results [6] and, in particular, the tightest 
constraint [7,8] 

-3.7 x l(r 17 /yr < d(\na)/dt < 3.1 x l(T 17 /yr (2) 

rely on the assumption that during these 2 billion years the value of a changed uniformly, if changed 
at all. This assumption looks rather natural but actually follows from nowhere. 

Thus in the modern epoch, at least on Earth since the Oklo times, the parameter a did not 
change more rapidly than by approximately 10~ 17 per year. If, on the other hand, we use the 
distant quasar data and take a mean value of d(lna)/dt for about 10 billion years, we shall obtain 
a variation rate of about 10~ 15 per year. Therefore one can conclude that at times earlier (maybe 
much earlier) than 2 billion years ago the value of a changed relatively rapidly but afterwards 
stopped or almost stopped to change. The task of theory was to explain such a behavior; however, 
if one takes into account the most recent observations [2,9], one should add the necessity of exlaining 
the spatial variations of a. Though, one cannot exclude the opportunity that the variations of a 
are purely spatial in nature whereas the time dependence is related to the fmiteness of the velocity 
of light: being located at a fixed point and at fixed time, we receive signals from distant regions 
of the Universe emitted at earlier cosmological epochs, and it is therefore impossible to separate 
spatial and temporal dependences of the parameters. 

Let us briefly discuss the theoretical models describing variations of a. Thus, following the 
pioneering ideas of Dirac and Eddington, Dicke and Peebles [10] in 1962 considered variations of 
a in cosmological models admitting a variable gravitational interaction intensity. Staniukovich [11] 
in 1965 discussed different variants of combined FPC variations in connection with Dirac's Large 
Number Hypothesis. Bekenstein [12] in 1982 described a model of a variations on the basis of the 
most general assumptions on the electromagnetic interaction: covariance, gauge invariance, causality 
and invariance with respect to time reversion. This led to a modified Maxwell electrodynamics and 
provided a certain dynamics of a. 

Since the advent of astronomical evidence on possible time variations of a, there emerged a 
whole class of new models describing such variations by introducing certain scalar fields. Thus, 
Sandvik et al. [13] proposed a cosmological extension of Bekenstein's theory [12] with a term of 
the form -iF^F^e" 2 ^ in the initial La grangian, where the scalar field ip interacts only with the 
electromagnetic field F^ u . The effect of the field ip in the dynamics of the expanding Universe 
was also considered. It was shown that in this model a remained almost constant in the radiation- 
dominated epoch, slightly increased in the matter-dominated epoch and approaches a constant value 
at times when the Universe expansion accelerates due to the presence of a positive cosmological 
constant. 

Spatial variations of a were also discussed much before they were claimed to be really discovered, 
in attempts to explain the discrepancy between cosmological and terrestrial data on a(t) [14,15]. 

In the recent attempts to explain both temporal and spatial variations of the fine structure 
constant, quite popular are models assuming the existence of domain walls connected with scalar 
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field dynamics (see, e.g., [16,17]). Thus, in [17] the initial action contains a dilaton-like scalar field 
interacting with the electromagnetic field and having a potential of the form V(<f>) = \\((fr 2 — rj 1 ) 2 . 
A domain wall is formed due to spontaneous symmetry breakdown. At points separated by the 
domain wall the values of a are different, which can explaing the observable variations if the wall 
intersects our Hubble volume. 

In [18] it has been shown that in F(R) gravity it is possible to obtain a static solution in the 
form of an effective (gravitational) domain wall, and that the choice of a logarithmic nonminimal 
interaction of the electromagnetic field with gravity in the form 



1 + ln 



R_ 

(where R$ is the modern value of the scalar curvature) makes it possible to describe variations of 
a, whose value grows as the curvature R decreases. 

Olive et al. [19] discuss a model with two domain walls, where the scalar field potential has three 
minima: 



V(<j>) = A 



(\®\ 2 - - V2ie ($ 3 - ($*) 3 ) + V 



It turns out that such a model much better describes the observational data than a similar one [17] 
with a single domain wall. 

The paper [20] suggests an extension of the previous BSBM (Bekenstein-Sandvik-Barrow- 
Magueijo) theory ( [12,13]) by introducing a dependence of the coupling constant u of the scalar 
field if) on the field itself, so that the Lagrangian contains the terms = —^(^d^ipd^ip and 
L cm = — ^F^F^e -2 ^. The choice of uj(i/j) allows for obtaining both growing and falling time 
dependences of a . This model differs from those with domain walls in that the variations of a are 
smooth and continuous, and a choice between these models must be easy with future more precise 
and reliable observational data. 

Mariano and Perivolaropoulos [21] have reported on a correlation between the spatial distri- 
bution of a values and the dipole anisotropy of the dark energy distribution. In the same paper 
they have suggested a theoretical model explaining this correlation (named "extended topological 
quintessence" ) which naturally predicts inhomogeneous spherical distributions of both the dark en- 
ergy and the values of a . The model assumes the existence of a huge global monopole with a size of 
Hubble order, which nonminimally interacts with the electromagnetic field. There emerge mutually 
related distributions of different parameters with a dipole anisotropy from the viewpoint of any 
observer located outside the monopole center. The monopole is formed after a phase transition in 
a set of three scalar fields with an 0(3) symmetric Lagrangian. 

In a later paper [22] the same authors support their inferences by the data on one more anisotropy 
also seeming to exist and to be aligned with other "dipoles", the so-called Large-Scale Velocity 
Flows (Dark Flow), i.e., recent indications that there is a large-scale peculiar velocity flow with an 
amplitude larger than 400 km/s on scales up to lOO/i" 1 Mpc (z < 0.03). 

It is also important to mention the theoretical models considering FPC variations in the frame- 
work of unification scenarios. In particular, P. Langacker et al. [23] consider possibile variations of 
coupling constants due to physics at very high energies, where the gauge couplings may be unified. 
It means that one should treat a joint variation of the fine structure and strong coupling constants. 

Similarly to [23], X. Calmet and H. Fritzsch [24] discuss FPC variations in the context of Grand 
unification. They show that such a consideration leads to small time shifts of the nucleon mass, 
the magnetic moment of the nucleon and the weak coupling constant, and it is expected to have a 
relative change of the nucleon mass larger than that of a by a factor of ~ 40 . 
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A more detailed discussion of theoretical models involving unification scenarios to explain the 
FPC variations can be found in [25]. 

It should be noted that all the above approaches, to explain variations of a, introduce scalar 
fields whose existence and manner of interaction with the electromagnetic field are postulated from 
the outset and are not explained in any way. In what follows, it will be shown how the scalar fields 
and their interaction law with electromagnetism naturally follow from nonlinear multidimensional 
gravity. Spatial variations of a are explained by a large-scale inhomogeneity of this scalar field. 
The magnitude of this inhomogeneity is constrained by CMB observations [26]. " 

The approach we are using has been formulated in [27], where a methodology was suggested 
allowing for a transition from multidimensional gravity with higher derivatives to Einstein-Hilbert 
gravity with effective scalar fields. Later on this approach was successfully applied for a unified 
description of the inflationary stage of the Universe and the modern secondary inflation [28] and an 
explanation of the origin of the Higgs field [29]; a mechanism of cascade reduction of multidimen- 
sional space to the observable one was suggested [30,31]. It has been shown under which conditions 
the compact extra dimensions become stationary (i.e., have a constant volume), and the cause of 
their maximum symmetry was found [32]. 

The present study has been performed in the framework of this approach and is an example of 
its employment. The paper is organized as follows. Sec. 2 briefly describes the general formalism 
used. In this framework, in Sec. 3 we build a homogeneous and isotropic cosmological model able 
to describe the present accelerated Universe along with a time dependence of the fine structure 
constant a. In Sec. 4, this cosmological model is slightly perturbed on large scale, which enables us 
to explain spatial variations of a. Sec. 5 is a brief conclusion. 



2 Multidimensional gravity and its reduction 

Consider a (D — 4 + di) -dimensional manifold with the metric 

ds 2 = g^dx^dx" + e 2f3{x) b ab dx a dx b (3) 

where the extra-dimensional metric components b a b are independent of x M , the observable four 
space-time coordinates. 

The D -dimensional Riemann tensor has the nonzero components 

BT = T? v 

± L pa A L pa 5 

R ab cd = e- 2 ^R ab cd + 5 ab cd ^, (4) 



where capital Latin indices cover all D coordinates, the bar marks quantities obtained from g^ v 
and b a b taken separately, (3^ = d^(3 and 5 ab c d = 5^8 b d — 5 d 5 b . The nonzero components of the Ricci 
tensor and the scalar curvature are 

Rp = Rp + di B^, 

R b a = e- 2 m b a + 6 b a [D/3 + d 1 (d^n 

R = R[ g ] + e" 2/r R[&] + 2d x up + dx{dx + l)(<9/3) 2 , (5) 

where (df3) 2 = (3^(3^ , □ = V M V M is the d'Alembert operator while R[g] and R[b] are the Ricci 
scalars corresponding to g^ and b a b, respectively. Let us also present, using similar notations, the 
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expressions for two more curvature invariants, the Ricci tensor squared and the Kretschmann scalar 
/C = R ABCD Rabcd '■ 

RabR AB = RfiuR 11 + 2d\R llv B^ u + d\B tlv B iJ ' h ' + e~^RabR 

+ 2e- 2/r R[6][n/3 + d l (d{5) 2 } + + d^dfi) 2 } 2 , (6) 

/C = + AdiB^B^ + e~^lC[b} +4e- 2 PR[b]{d(3) 2 + 2d 1 {d x - l)[(d/3) 2 ] 2 . (7) 

Suppose now that b ab describes a compact d\ -dimensional space of nonzero constant curvature, 
i.e., a sphere (K = 1) or a compact d\ -dimensional hyperbolic space (K = —1) with a fixed 
curvature radius r normalized to the D -dimensional analogue m D of the Planck mass, i.e., r = 
l/m D (we use the natural units, with the speed of light c and Planck's constant h equal to unity). 
We have 

R a \ d = Km 2 D 5 ab cd , 
R b a = Km 2 D {d l -l)5l 

R[b] = Km 2 D d x {d x - 1) = R b . (8) 

The scale factor b(x) = e@ in (3) is thus kept dimensionless; Rb has the meaning of a characteristic 
curvature scale of the extra dimensions. 

Consider, in the above geometry, a sufficiently general curvature-nonlinear theory of gravity with 
the action 



S = l m °~ 2 J V°9d D x(L g + L m ), 



L g = F(R) + Cl R AB R AB + c 2 /C, (9) 

where F(R) is an arbitrary smooth function, c\ and c 2 are constants, L m is a matter Lagrangian 
and D g = \ det(g MN ) \ . 

The extra coordinates are easily integrated out, reducing the action to four dimensions: 

S = l -V[d l ]m 2 D j ^d A xe d ^[L g + L m ], (10) 

where A g = \ det((7 AiI/ )| and V[d\] is the volume of a compact d\ -dimensional space of unit curvature. 

Eq. (10) describes a curvature-nonlinear theory with non-minimal coupling between the effective 
scalar field /3 and the curvature. Let us simplify it in the following way (putting, for convenience, 
m D = 1, so that all quantities are now expressed in (D -dimensional) Planck units: 

(a) Express everything in terms of 4D variables and (3(x); we have, in particular, 

R = R4 + (f) + f u R, = R[g], f 1 = 2d 1 nf3 + d 1 (d 1 + l)(d(3) 2 , (11) 
where we have introduced the effective scalar field 

<f>(x) = R b e~ mx) = Kd^dx - 1) e" 2/3W (12) 

The sign of <ft coincides with k — ±1, the sign of curvature in the d\ extra dimensions. 

(b) Suppose that all quantities are slowly varying, i.e., consider each derivative <9 M (including those 
in the definition of R) as an expression containing a small parameter e; neglect all quantities of 
orders higher than 0(e 2 ) (see [27,31]). 
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(c) Perform a conformal mapping leading to the Einstein conformal frame, where the 4-curvature 
appears to be minimally coupled to the scalar 0. 

In the decomposition (11), both terms f\ and R 4 are regarded small in our approach, which 
actually means that all quantities, including the 4D curvature, are small as compared with the 
D -dimensional Planck scale. The only term which is not small is 0, and we can use a Taylor 
decomposition of the function F(R) = F(<f> + R 4 + fi) : 

F(R) = F(4> + R 4 + /i) ~ F(cf>) + F'(cf>) ■ (R 4 + }\) + (13) 

with F'((p) = dF/d(f>. Substituting this, and the corresponding decompositions of the expressions 
(6) and (7), into Eq. (10), we obtain, up to 0(e 2 ), the following effective gravitational Lagrangian 
L g in Eq. (10): 

L g = F'(4>)R 4 + F{<j>) + F'{<P)h + c*4> 2 + 2ci0D/3 + 2{cyd x + 2c 2 )(<9/3) 2 (14) 

with c* = Cx/di + 2c 2 /[d 1 (di - 1)]. 

The action (10) with (14) is typical of a scalar-tensor theory (STT) of gravity in a Jordan frame. 
To study the dynamics of the system, it is helpful to pass on to the Einstein frame. Applying the 
conformal mapping 

g,u ^ g,u = \f(<i>)\g^ f{<f>) = e dl/3 F'(0), (15) 

after a lengthy calculation, we obtain the action in the Einstein frame as 



S= iv[di] J v^(signF')L 



L = R 4 + Kv{<p){d<PY - 2^e(0) + L m , (16) 



e -diP 



L m = {signF') Ji ^ 2 L m] (17) 
'F"\ 2 , ^ , 4(c 1 + c 2 ), 



F' 



e -di/3 



-2VM = (signF')— p[F(0) +c*0 2 ], (19) 

where the tilde marks quantities obtained from or with 7j^ v ; the indices are raised and lowered with 
g^,; everywhere F = F(4>) and F' = dF/d<j); e@ is expressed in terms of using (12). 

Let us consider the electromagnetic field F^ v as matter in the initial Lagrangian, putting 

L m = a^F^F^, (20) 

where ot\ is a constant. After reduction to four dimensions this expression acquires the factor e dl/3 
arising from the metric determinant: \[~®g~ = y^ge dl/3 . In the subsequent transition to the Einstein 
picture the expression \f*gF tlv F tw remains the same (it is the well-known conformal invariance of 
the electromagnetic field), hence the Lagrangian (17) takes the form 

L m = a^e d ^F^, (21) 

and for the effective fine structure constant a we obtain 

— = e dl ^-«, (22) 

where «o and (3 are values of the respective quantities at a fixed space-time point, for instance, 
where and when the observation is taking place. 
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3 The cosmological model 

Depending on the choice of F(R) , the parameter c\ and c 2 and the matter Lagrangian in the action 
(9), the theory under consideration can lead to a great variety of cosmological models. Some of 
them were discussed in [27], mostly those related to minima of the effective potential (19) at nonzero 
values of 0. Such minima correspond to stationary states of the scalar and consequently of the 
scale factor b = e^ of the extra dimensions. If the minimum value of the potential is positive, it 
can play the role of a cosmological constant that launches an accelerated expansion of the Universe. 

Here, we would like to focus on another minimum of the potential Vain, existing for generic 
choices of the function F(R) with F' > and located at the point = 0. The asymptotic — > 
corresponds to growing rather than stabilized extra dimensions: b = e@ ~ 1/ y|0[ — > oo. A model 
with such an asymptotic growth at late times may still be of interest if the growth is sufficiently 
slow and the size b does not reach detectable values by now. Let us recall that the admissible range 
of such growth comprises as many as 16 orders of magnitudes if the D -dimensional Planck length 
l/m D coincides with the 4D one, i.e., about 10 -33 cm: the upper bound corresponds to lengths 
about 10 -17 cm or energies of the order of a few TeV. This estimate certainly changes if there is no 
such coincidence. 

One should note that small values of to be considered here are still very large as compared to 
4D quantities, and so our general assumptions are well justified. Indeed, according to (12), 

di(di — 1) 



b 2 > 

where b < 10 16 , hence |0| > d\ ■ 10~ 32 , while the quantity R4, if identified with the curvature of 
the modern Universe, is of the order 10~ 122 in Planck units (that is, close to the Hubble parameter 
squared, or (the Hubble time) -2 , see also Eq. (34) below). 

Let us check whether it is possible to describe the modern state of the Universe by an asymptotic 
form of the solution for — > as a spatially flat cosmology with the 4D Einstein-frame metric 

d^\ = dt 2 -a 2 (t)dx 2 , (23) 

where a(t) is the Einstein-frame scale factor. We shall be are working in the framework of quadratic 
gravity with a cosmological constant, i.e., 5 

F(0) = -2A D + F 2 2 , (24) 

where Ad is the initial cosmological constant. Then, substituting F' = 20 and F" = 2, we obtain 
for the kinetic and potential terms in the Lagrangian (16) in the first approximation in 0: 



K E « K /(2(f> 2 ), K = - -df-d 1 + Q + 2(ci + c 2 



1 



2 

A, 



V^V e-^, Vo-^f^, 23=* -4. (25) 

It is clear that this model can work only if d\ > 4. In terms of /3 instead of 0, the Lagrangian 
takes the form 

L = R 4 + 2K (df3) 2 - 2V e~ 2 ~ dp + L m , (26) 



5 Wc assume for certainty </> > 0, or, which is the same according to (12), K = +1, but everything can be easily 
reformulated for d> < 0. 
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Neglecting the gravitational influence of the electromagnetic field (that is, considering only vacuum 
models), one can write down the independent components of the Einstein and scalar field equations 
with the unknowns (3(t) and a(t) in the form 



a 2 



3^ = K (3' 2 + V e- 2d ^ (27) 



a 2 

H + 3 ^ = Y^L e -^. (28) 
a K 

These equations, corresponding to a scalar field with an exponential potential, can be solved 
exactly but the solution looks rather involved, and for our purpose more preferable is the compar- 
atively simple approximate solution that can be obtained in the slow-rolling approximation; the 
latter should be acceptable at late times. Let us suppose that 

|/3|«3^/3, K /3 2 «K.e- 2 ^, (29) 

and neglect the corresponding terms in Eqs. (27) and (28). Then, expressing the quantity a/ a from 
(27) and substituting it into (28), we obtain 

/3 = ^e-^, (30) 



whence 



-2 



d IVo, 



e '=^Vf (i+ii) - (31) 

where t\ is an integration constant. For the scale factor a(t) we have 

— = — - — =>■ a = a 1 (t + ti) p , ai = const, P—^- (32) 
a t + ti d 

Substituting the solution to the slow- rolling conditions (29), we make sure that they hold as 
long as p 3> 1 , or in terms of the input parameters of the theory, 

d 2 - 2d 1 + 12 + 4( Cl + c 2 ) 

v = (d^W y> L (33) 

We will assume that this condition holds. 

A further interpretation of the results depends on which conformal frame is regarded physical 
(observational) [33,34], and this in turn depends on the manner in which fermions appear in the 
(so far unknown) underlying unification theory involving all interactions. 

Let us adopt the simplest hypothesis that the observational picture coincides with the Einstein 
picture and make some estimates. Thus, the inverse of the modern value of the Hubble parameter 
(the Hubble time) is estimated as 

t H = 1/H = a /d w 4, 4 x 10 17 c w 8 x 10 60 t ph (34) 

where t p \ is the Planck time and the index "0" marks quantities belonging to the present time, 
which is a usual notation in cosmology. From (32) it follows that H = p/(to + ti), whence 



U ■= t + h = pt H » t H - 



(35) 
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With p » 1, the model satisfies the observational constraints on the factor w in the effective 
equation of state p = wp of dark energy that causes the accelerated expansion of the Universe: at 
w = const we have a ~ ^2/(3+3^) ^ consequently, w = —1 + 2/(3p) is a number close to —1: for 
example, to have w ~ —0.99, one should put only p = 66. Meanwhile, the recent observational data 
allow for a comparatively large range of w [35-38] but anyway admitting w = —1 corresponding to 
a cosmological constant. This follows from combining the recent measurements of cosmic microwave 
background anisotropies, Supernovae luminosity distances, baryonic acoustic oscillations, and H (z) 
measurements, though different tests lead to different confidence intervals. 

Furthermore, the "internal" scale factor b(t) = e^ grows much slower than a(t): 

m-^f\ «) 

Using the expression for Vq from (25), one can estimate the initial parameter A^, connecting it 
with the size of the extra factor space 60 : i n Planck units, 

A D = \2HldKd, - « ^d\{d x - l) 2 ^" 4 x 1(T 120 . (37) 

16 

As already mentioned, the "internal" scale factor b = e 13 should be in the range 1 <^ bo < 10 16 
in Planck units. The estimate (37) shows that the present model makes much easier the well- 
known "cosmological constant problem" (the difficulty of explaining why in standard cosmology 
^-standard ~ 1CT 122 in Planck units). For instance, if (in the admissible range) b = 10 15 and 
di — 12, it follows An = 3267 without any indication of fine tuning. 

Let us estimate the possible range of the parameters C\ and ci in the action (9). The present 
model describes only the modern stage of the Universe evolution, but it should admit an improve- 
ment after which it will account for other stages, including the early inflation. Then one should 
require that the curvature-nonlinear terms in the initial Lagrangian should not violate our slow- 
change approximation, see Sec. 2 This leads to the condition <C 10 11 . Indeed, during inflation, 
the Hubble parameter is H ~ 10~ 6 in Planck units, while the scalar curvature at inflation, when 
the 4D geometry is approximately de Sitter, is estimated as R ~ 12H 2 ~ 10 -11 . Assuming that the 
Ricci and Riemann tensor components have the same order of magnitude, we find that the condition 
R ^> C\R AB Rab , used above in the framework of the slow-change approximation, will be violated 
if Ci is too large. The upper bound of the parameter C2 is obtained in a similar way. 

The smallness of the observed variations of a leads to another constraint on c\ and C2 : according 
to (38), 

/ (hlh\-d* / t + ti N- 2 W-4) 2d 1 t-t 

so that a/a ~ 10 _10 /p per year. By the empirical data, this quantity cannot be larger than about 
10~ 17 per year. A comparison leads to the constraint p > 10 7 and hence the effective equation-of- 
state parameter w is equal to —1 up to seven meaningful digits. Taking into account the relation 
(33) between p and the input parameters c\ and C2 , we obtain similar bounds on these parameters 
if the number of extra dimensions d\ is not too large. 

Thus the allowed range of C\ and ci (assuming that they are of the same order of magnitude), 

10 7 < ci, 2 < 10 11 (39) 



is wide enough, which means that any fine tuning is absent. 
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One of the well-known constraints on Kaluza-Klein-like cosmologies is the requirement of a 
sufficiently slow evolution of the internal scale factor, e.g., according to [39], the extra- dimensional 
volume should not have changed by more than 10 % since the times of primordial nucleosynthesis. 
This requirement rests on the relation between the effective Newtonian gravitational constant 
and the volume of extra dimensions. However, in the Einstein conformal frame used here, = 
const by definition, therefore the above constraint does not apply here. Still, even if we passed on 
to the Jordan frame, the relative variations of G^ would be the same as those of a and obey the 
law (38), i.e., within about 10~ 5 for the Hubble time. 

In the next section we shall see that the inequality p > 10 7 and consequently > 10 7 are 
substantially relaxed in the perturbed model. 



4 Spatial variations of a 

In the previous section we discussed the properties of a homogeneous model which does not contain 
any spatial variation of a (and any other physical quantity). Let us try to describe variations of 
a by taking into account spatial perturbations of the scalar field and the metric. Only long-wave 
perturbations will be of interest for us, with characteristic lengths of the order of the horizon size. 

An observed statistically isotropic sky means that there is no preferred axis. Nevertheless, 
super-horizon components were produced by quantum fluctuations at the beginning of inflation in 
the same way as fluctuations of smaller scale. It means that the dipole component must exist though 
hardly observed due to its contamination by the Doppler effect caused by the motion of our Local 
Group with respect to the CMB. Hence there must exist a weakly expressed distinguished direction 
along which the metric and scalar field inhomogeneity is most clearly pronounced. 

Accordingly, we now choose a metric more general than (23), 

dsl = e 2S ^dt 2 - a(t) 2 e 25x dx 2 - a(t) 2 e 25 \dy 2 + dz 2 ), (40) 

where x is the distinguished direction and £7, 5\, 5t] 1 are functions of x and t. In addition, we 
replace the effective scalar field (3{t) with (3(t) + 5(3(x,t). 

Then the relevant Einstein-scalar equations corresponding to the Lagrangian (26) can be written 
as follows (preserving only terms linear in the "deltas"): 

5(3 + -8$ + $(6X - &y) - - 2 5f3" + ^-5(V? e 2 ^) = 0, (41) 

a (5A-5 7 ) = 5(Ve 2 "'), (42) 

5 1 ' = K f3 5(3', (43) 



a 
a 

a 



where we have chosen the gauge (in other words, the reference frame in perturbed space-time) 
5r] = 0, the dot and the prime stand for d/dt and d/dx, respectively. We have also denoted 
V = V E = V e~ 2l(3 and V p = dV/d/3. 

Integration of (43), without loss of generality, leads to 

5 7 = (44) 

where, as before, H = a/ a. This equation enables us to estimate the quantity 5/3 . Indeed, according 
to the CMB data [26], we can take £7 ~ 10~ 5 , while the coefficient before 5(3 is of the order of 
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unity at the present epoch (according to Eq. (32) we have (K /H )f3(t = t ) = K /(d-p) = d ~ 1), 
we obtain 5(3 ~ 10~ 5 . 

Substituting this £7 to (41) and taking the difference 5X — £7 from (42), we finally arrive at the 
following single wave equation for 5(3: 



5(3 + —5P - -M" + 5(3 
a er 



0. (45) 



with an arbitrary constant Kq and an arbitrary potential V(/3). In our case, with V = V^e 2d/3 
and Kq given in (25), we obtain 

6P + ™ Sj 3 - ^5(3" + ^ 5(3 = 0, (46) 

a a 1 p 

while the background quantities a(t) and (3{t) are determined by the solution (31), (32). It remains 
to find a solution for 5(3 which, being added to the background (3(t), would be able to account for 
the observed picture of variations of a. 

Since the background is ^-independent, we can separate the variables and assume 

5(3 = y(t) sin k(x + xo) 

where k has the meaning of a wave number, of order of the cosmo logical horizon scale, and y(t) 
must be as small as 10~ 5 . Then y(t) obeys the equation 



3p . 

y + TTu y + 



k 2 Qp 
+ 



ai(t + h) 2 P (t + ti) 5 



y = 0. (47) 



Since the equation (47) has been derived in a certain approximation and describes only a restricted 
period of time close to the present epoch, it is reasonable to seek the solution in the form of a Taylor 
series: 

y(t) = y + yi(t -t ) + ^y 2 (t -t ) 2 + . . . , yi = const. (48) 

Then yo and y± can be fixed at will as initial conditions, and Eq. (47) leads to expressions of 
2/2, Hz-, •• • m terms of yo and y±. Even more than that, for a certain neighborhood of t — to 
we can simply suppose y = y + yi(t — t ) ■ Actually, this approximation is good enough for 
t — to <C = to — ^1 • 

In this approximation we obtain the following expression for variations of a: 

H « 1 _ ^ tzh - dl S in[k(x + xo)] [y + Vi(t - t )} + 0(e 2 ), (49) 

where 0(e 2 ) means 0((t — to) 2 /t 2 ). Assuming that the observer is located at x = and requiring 
a/a = 1 + 0(e 2 ) at x = 0, we obtain the condition 

y 1 sm(kxo) = —l/(dU). (50) 

This explains very small, if any, variations of a on Earth at present and since the Oklo times. 
Indeed, since 2xl0 9 yr »s ^tn while = ptn, the addition 0(e 2 ) is of the order of l/(50p 2 ), where 
p > 1. If we take, for instance, p = 1000, then at the Oklo time (2xl0 9 years ago) we obtain a 
relative a variation of the order 0.5xl0~ 8 , which makes about 0.25 xlO -17 per year. 
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A substitution of (50) and (50) into (49) at t — to = —x for x > gives 

a/ao ~ 1 — diyo sm(kxo) + diy kx cos(/cxo) + 0(e 2 ) (51) 

at x <C t* . The same result is obtained if we substitute t — t = x for x < . 

Fig. 1 compares the observational data and the predictions of our model with the parameters 
indicated there. Recall that we are considering long-wave fluctuations, such that k < 1/rjj ~ 0.1 
(billion years)" 1 (where r# is the modern horizon size), with a small magnitude yo < 10~ 5 . The 
relations obtained are in good agreement with these estimates. We are using the conventional 
normalization ao = 1 . 




Figure 1 : The r dependence of 5a / ao (the distance r is measured in billions of light years) . The 
dashed lines correspond to Eq. (1), the solid red line to Eq. (49) at the parameter values d\ = 12, 
p = 10 7 , y = -4.7 xlO" 6 , yt = -10" 7 (bill, years)" 1 , k = 0.02 (bill, light years)" 1 , x = 1 
billion of light years. 

Evidently, our model, in addition to the input theoretical parameters like di, Ci, c 2 , contains 
the parameters k, Xq, yo, y±, depending on the initial form of the extra space metric. 

In the framework of chaotic inflation, these parameters vary in different regions of the visible 
part of the Universe. Their choice enables us to explain the spatial variations of a in agreement 
with the observations [2]. Actually, there are only two conditions imposed on them: (50) and the 
relationship identifying (51) with the expression (1) at r = x and cosip = 1, i.e., on the dipole axis. 
We obtain (in Planck units) 

dmk cos(kx ) w -2 x 10" 66 . (52) 

(This numerical value is used for obtaining the solid line in Fig. 1.) The small constant shift of 
the a value at x = against the background does not change the interpretation of the results 
obtained. It should be stressed that (52) is not a fine-tuning relation but simply fitting of the model 
parameters to the observational data. In fact, the very small number in the r.h.s. of (52) results 
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from the natural scale of k ~ 1/rzj, where r# ~ 10 28 cm ~ 10 61 £ p \ is the Hubble radius in terms 
of the Planck length £ p \; five more orders of magnitude in (52) are related to the smallness of a 
variations. 

The input parameters C\ and c 2 are now not so strongly constrained by the condition of slow 
variations of a on Earth: this condition is already provided by the equality (50) if we take p > 1000, 
hence C\ + c 2 > 1000. The approximation p ^> 1, in which our solution has been obtained, then 
also holds quite well. The inequality (39) is thus replaced by a much weaker one: 

1000 < cj ~ c 2 < 10 11 . (53) 

5 Conclusion 

We have studied the possible effect of extra dimensions on large-scale variations of the fine structure 
constant a in space and time. In the multidimensional paradigm under consideration, the observable 
values of a and probably other physical quantities, including fundamental constants, depend on 
the size of the extra factor space. Variations of the dark energy density can be mentioned as an 
example. Indeed, the space-time variations of the energy density are dominated by those of the 
potential V = Ve given in (25). The relative variation 5V/V = —2d5f3 is of the same order of 
magnitude as the space-time variations of a according to (38). They are too small to be observed 
in the near future. 

We have discussed the dipole component only, but it seems evident that the same basis is 
applicable to higher multipoles in a variations. It means that the observational data, being quite 
uncertain, "feel" these components, and further observations may detect them. 

We have focused on the behavior of a because it is the only fundamental constant for which 
there are more or less reliable data indicating its variations. We are also planning to analyze the 
behavior of other constants, above all, the gravitational constant and the particle masses. 

The model described here does not consistently include other kinds of matter than dark energy 
(represented by a scalar field of multidimensional origin). However, even such a simple model shows 
an agreement with the observational data (see Fig. 1). The same numerical parameters also well 
agree with the CMB constraints which impose an upper bound on the fluctuation magnitude of the 
extra-dimensional metric. 

An advantage of the present model of a variation against many others (e.g., [16,17,19]) is that it 
assumes a common origin of dark energy and FPC variations. Our model also predicts the existence 
of higher multipoles in a variations. 
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